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ABSTRACT 

We derive a variety of physical parameters including star formation rates (SFRs), dust attenuation and burst 
mass fractions for 6472 galaxies observed by the Galaxy Evolution Explorer ( GALEX) and present in the SDSS 
DR1 main spectroscopic sample. Parameters are estimated in a statistical way by comparing each observed 
broad-band SED (two GALEX and five SDSS bands) with an extensive library of model galaxy SEDs, which 
cover a wide range of star formation histories and include stochastic starbursts. We compare the constraints 
derived using SDSS bands only with those derived using the combination of SDSS and GALEX photometry. We 
find that the addition of the GALEX bands leads to significant improvement in the estimation of both the dust 
optical depth and the star formation rate over timescales of 100 Myr to 1 Gyr in a galaxy. We are sensitive to 
SFRs as low as 1 0~ 3 M© yr" 1 , and we find that low levels of star formation (SF) are mostly associated with early- 
type, red galaxies. The least massive galaxies have ratios of current to past-averaged SF rates (^-parameter) 
consistent with constant SF over a Hubble time. For late-type galaxies, this ratio on average decreases with 
mass. We find that b correlates tightly with NUV -r color, implying that the SF history of a galaxy can be 
constrained on the basis of the NUV -r color alone. The fraction of galaxies that have undergone a significant 
starburst episode within the last 1 Gyr steeply declines with mass — from ~ 20% for galaxies with ~ 10 8 M Q to 
~ 5% for ~ 10 u Mq galaxies. 

Subject headings: galaxies: evolution- galaxies: fundamental parameters — galaxies: starburst — ultraviolet: 
galaxies 



1. INTRODUCTION 

Modern large-scale galaxy surveys are allowing new con- 
straints to be placed on the history of star formation (SF) in 
galaxies. High-quality optical spectra collected by the Sloan 
Digital Sky Survey (SDSS) have been used to study the re- 
cent star formation histories, dus t content and m etallicities 
of over 10 s nearby galaxies (e.g., Kauffman n et alJ l2003a b; 
iBrinchmann et alJl2004t iTremonti et al.ll2004l) . These analy- 
ses make use of specific absorption and emission lines in the 
galaxy spectra and employ new models of the spectral evolu- 
tion of galaxies, which include a physically consistent treat- 
ment of the production of starlight and it s transfer through the 
interstellar medium (ISM) in galaxies ( Chariot & Falll l2000t 
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ICharlot & Longhettill200ll iBruzual & Chario3l200l . In ad- 
dition, the modeling accounts for the stochastic nature of SF 
when interpreting galaxy spectra by using large Monte Carlo 
libraries of different SF histories to estimate physical param- 
eters such as stellar mass, age, and SF rates in a statistical 
fashion. 

In this Letter, we use a similar approach to interpret the 
combined ultraviolet (UV) and optical colors of 6472 SDSS 
galaxi es observed b y the Galaxy Evolution Explorer ( GALEX) 
JMartin et aTll2004l) . We show that the addition of UV infor- 
mation to the optical SEDs of galaxies leads to significant im- 
provements in the estimates of the star formation rates (SFRs), 
starburst histories and dust attenuations. 

2. DATA AND SAMPLE 

2.1. GALEX and SDSS data 

We consider galaxies with combined GALEX and SDSS 
photometry, for which spectroscopic redshifts are avail- 
able from the SDSS. GALEX images the sky at far-UV 
(FUV; 1530 A) and near-UV (NUV; 2310 A) bands in two 
modes: All-sky Imaging Survey (AIS, raii m (AB) m 20.4) and 
Medium-deep Imaging Survey (MIS, raii m (AB) « 22.7). Each 
circular field covers 1.1 sq. deg. Here we use the internal re- 
lease of the catalog (IR0.2) that consists of 649 AIS and 94 
MIS fields, of which 117 and 91 overl ap (at least in part) with 
the SDSS Data Release One (DR1, lAbazaiian etal] 120031) . 
We use FUV and NUV fluxes and errors derived in elliptical 
apertures. For GALEX sources with counterparts in the SDSS 
DR1 spectroscopic sample we extract SDSS ugriz model col- 
ors normalized to the Petrosian r magnitude, 90% and 50% 
Petrosian /-band radii, and the likelihoods of exponential and 
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TABLE 1 

Average Galaxy Parameters and Average Parameter Errors 



Parameter name 


(G + S) 


RMS A(G + S) a 


<<t(G + S)) 


(o-(S)) 


Gain 


t i 
'gal 


6.50 


0.53 


1.91 


2.16 


12% 


7 


0.40 


0.06 


0.21 


0.22 


7% 


Z/Zq 


1.00 


0.14 


0.36 


0.42 


12% 


TV 


1.24 


0.17 


0.60 


0.86 


29% 


I" 


0.36 


0.03 


0.14 


0.15 


6% 


A-FUV 


2.03 


0.23 


0.62 


1.05 


41% 


Anuv 


1.46 


0.18 


0.46 


0.78 


41% 


logM, 


10.42 


0.02 


0.08 


0.09 


11% 


log SFR(100 Myr) 


0.00 


0.07 


0.27 


0.54 


49% 


log SFR(1 Gyr) 


0.18 


0.12 


0.30 


0.55 


45% 


^burst 


2.31 


0.90 


1.62 


1.65 


2% 


F burst (100 Myr) 


0.00 


0.06 c 


0.01 


0.01 


39% 


fburstO Gyr) 


0.02 


0.08 c 


0.04 


0.06 


39 c /r 
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FIG. 1 . — Comparison of galaxy stellar masses derived in Kauffmann et al. 
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NOTE. — f ga i , /burst and 7 1 are in Gyr, SFR in M@ yr 1 , M, in Mq , and Af yv 



and A NUV in magnitudes. effective V-band absorption optical depth Ty affecting stars 

RMS change in the parameter value when changing the frequency of bursts from younger than 10 Myr (that arises from giant molecular clouds 

and the diffuse ISM) and the fraction fi of it contributed by 
the diffuse ISM, that also affects older stars. Based on pre- 



50% to 10% in the last 2 Gyr in the model libraries 
b Improvement in a by adding the GALEX fluxes to SDSS ugriz, constraints. 



c Only for galaxies for which Fburst Is > 0. 



de Vaucouleurs r-band light profiles. 

2.2. Sample 

We construct our sample from G ALEX objects match ed to 
SDSS DR1 objects in a 6" radius JSeibert et alJl2004l) . We 
then restrict the sample to objects spectroscopically classified 
as galaxies and with redshifts 0.005 < z < 0.25. To avoid 
mismatches or matching an unresolved GALEX object with 
resolved SDSS objects, we additionally impose astrometric 
criteria which remove rj 10% of the galaxies, leaving 6826. 
About 70% of the remaining galaxies have redshifts in a nar- 
row range between 0.05 and 0. 15 (the number of SDSS galax- 
ies with spectra drops steeply when z > 0.1). GALEX MIS 
and AIS recover on average 80% and 50% of the SDSS main 
galaxy spectroscopic sample {r\\ m < 17.8). MIS galaxies con- 
stitute 60% of our sample. The addition of UV selection leads 
to incompleteness for galaxies redder than u-r ps 3.0 in the 
MIS and u-r « 2.5 in the AIS. However, brighter red galaxies 
are still present in our sample. 

3. MODEL GALAXY PARAMETERS 

3.1. Libraries of galaxy models 

We use the lBruzual & Charlotl ( 120 03) population synthesis 
code to generate libraries of stochastic realizations of model 
SF histories at five redshifts equally spaced between 0.05 
and 0.25. Each l ibrary contains ~ 10 5 m odels. For sim- 
plicity, we follow Kauffmann et al. (2003a) and parameter- 
ize each SF history in terms of two components: an under- 
lying continuous model with an exponentially declining SF 
law (SFR[f] oc exp[~7?]) and random bursts superimposed on 
this continuous model. 12 Attenuati on by dust is descri bed by 
means of two main parameters (see Chariot & Fall 2000): the 

12 Specifically, we take the galaxy age f ea i to be distributed uniformly over 
the interval from zero to the age of the universe at a given redshift, and the SF 
timescale parameter over the interval < 7 < 1 Gyr~' . Random bursts occur 
with equal probability at all times until r ga | . They are parameterized in terms 
of the ratio between the mass of stars formed in the burst and the total mass 
of stars formed by the continuous model over the time f Ba i . The ratio is taken 
to be distributed logarithmically from 0.03 to 4.0. During a burst, stars form 
at a constant rate for a time distributed uniformly in the range 30-300 Myr. 



vious analyses of SPSS galaxies dKauffmann et alJ l2003a. 
iBrinchmann et alJl20 04). we take Ty to be distributed from 
to 6 with a broad peak around 1 and /j, to be distributed 
from 0. 1 to 1 with a broad peak around 0.3. Our model galax- 
ies have metallicities uniformly distributed between 0. 1 and 2 
Z Q . 

3.2. Deriving physical parameters 

We derive statistical estimates of various physical param- 
eters by comparing the observed SED of each galaxy in the 
sample to the model SEDs as follows. First, we convert SDSS 
magnitudes into fluxes, correcting for th e offsets between the 
SDSS- and AB-magnitude zeropoints ( Abaz aiian et al 120041) 
and applying small zeropoint offsets to GALEX fluxes. To 
SDSS flux errors we add calibration errors in quadrature. Af- 
ter accounting for Galactic reddening, we compare the SED of 
the observed galaxy to all SEDs in the closest-redshift model 
library. The \ 2 goodness of fit of each model determines the 
weight (oc exp[~x 2 /2]) to be assigned to the physical parame- 
ters of that model when building the probability distributions 
of the parameters of the given galaxy. The probability density 
function (PDF) of a given physical parameter is thus obtained 
from the distribution of the weights of all models in the li- 
brary. We characterize the PDF using the median and the 2.3- 
97.7 percentile range (equivalent to ±2er range for Gaussian 
distributions), and also record the y 2 of the best-fitting model. 
We perform this analysis for all galaxies in our sample. The 
distribution of \ 2 values of the best-fitting models is gener- 
ally very good, implying that our libraries do reproduce the 
observed SEDs. However, there is a tail of large \ 2 values- 
usually identifiable as objects with suspect SDSS magnitudes, 
or as "shredded" objects. We remove ~ 300 galaxies with 
poorest fits. In the end, we obtain estimates of physical pa- 
rameters for 6472 galaxies. 

3.3. Parameters and their errors 

We estimate physical parameters in two ways: using the 
GALEX FUV and NUV fluxes combined with the SDSS ugriz 
fluxes; and using only the SDSS fluxes. This allows us to 
quantify the effect of adding ultraviolet constraints. 

The burst probability is set so that 50% of the galaxies in the library have 
experienced a burst in the past 2 Gyr. 
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FIG. 2. — Dependence of metallicity, attenuation and star formation parameters on fundamental galaxy properties (see text). Blue dots have SDSS profiles 
that are more likely to be exponential (late-type galaxies), while red dots have higher likelihood of S 1 ' 4 profiles. (Note that this division may misclassify some 
galaxies.) 



Table 1 summarizes our results. For each physical param- 
eter, we list the sample-averaged parameter value ((G + S)) 
derived using GALEX + SDSS constraints, "la" estimate of 
the average parameter error with (G+S) and without (S) in- 
cluding the GALEX constraints, and the increase in accuracy 
(gain) achieved with GALEX . We also investigate the effect 
of changing the fraction (from 50% to 10%) of galaxies in the 
model libraries that had bursts over the past 2 Gyr. In Table 1, 
we report the resulting RMS change in parameter value. It is 
usually much smaller than the "fitting" error, except for the 
time fburst since the last burst of SF. 

The physical parameters for which the estimates are most 
significantly improved with GALEX are the 100 Myr- and 
1 Gyr-averaged SF rates, (SFR(100 Myr) and SFR(1 Gyr)), 
FUV and NUV dust attenuations (Afuv ar, d A-nuv), and es- 
timates of the fraction of a galaxy's stellar mass formed in 
bursts (Fbur S t(100 Myr) and fb urst (l Gyr)) over the last 0. 1 and 
1 Gyr (although it is sensitive to the assumed frequency of 
bursts in the model libraries). 

4. ANALYSIS OF PHYSICAL PROPERTIES, STAR FORMATION AND 
STARBURST HISTORIES 

In Figure ^ we compare the stellar masses derived from 
our analysis of GALEX + SDSS colors with those derived by 
IKauffmann et alj j2003a). These authors used the strengths of 
the HSa absorption-line index and 4000 A break in the SDSS 
spectra (that trace the recent and past-averaged SF histories) 
to improve the constraints on the stellar mass-to-light ratios 
relative to estimates based on a single optical color. The scat- 
ter in the difference between the two types of mass estimates 



(0.11 dex without 3a outliers) is very well matched by the 
uncerta inties in the two studies (0.08 dex). Kauffman n et alJ 
(2003a) also compute the fraction of a galaxy's stellar mass 
formed in bursts over the last 2 Gyr. In our study, bursts 
over 2 Gyr are poorly constrained. This is because the UV 
light is less sensitive to stars older than 1 Gyr than the HSa 
index. We have checked that galaxies for which we find 
^burst(l Gyr) > 0.05 have 400 A breaks and H5 A str engths 
typical of galaxies for which IKauffmann et alJ J2003al) find 
F burst (2 Gyr)> 0.05. 

We present our main results in Figure[2] where we show the 
dependence of selected physical parameters on three funda- 
mental galaxy properties: stellar mass, galaxy type, and color. 
We use the concentration parameter Rgo /R50 as a proxy for 
galaxy type. Values of R90/R50 larger than 2.5 correspond 
mostly to early-type galaxies, but late-ty pe galaxies can also 
have large Rgo/R^o (Fukugita et all 120041) . We choose the 
NUV - r color as it provides a broad baseline and is not sub- 
ject to large k-corrections. Objects are color-coded to indicate 
whether an exponential profile (blue) or a de Vaucouleurs 7? 1 / 4 
profile (red) has higher probability. This serves as an addi- 
tional rough indicator of galaxy type. 

The first row of panels in Figure [2] presents the derived 
stellar metallicities. It was recently demonstrated that star- 
forming galaxies in SDS S exhibit a tight m etallicity (of ISM) 
vs. mass correlation (Tremonti et al. 2004). Although metal- 
licity is not strongly constrained using our method (see Ta- 
ble 1), we do find that low-mass galaxies do not reach metal- 
licities as high as high-mass galaxies. (Note that galaxies with 
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FIG. 3. — Fraction of galaxies of a given mass with starbursts. Histograms 
correspond to fraction of galaxies with bursts (that produced more than 5% 
of stellar mass) having occurred in the last 100 Myr and 1 Gyr (dashed and 
filled lines) and . The fraction of galaxies with recent bursts declines sharply 
with galaxy mass. 



M* < 10 I0 M Q are preferentially disk galaxies, and mostly 
have NUV - r < 4.) Massive galaxies span a wide range in 
Z, with blue and red dots occupying the same space. 

The second row of panels shows that low-mass galaxies on 
average suffer less attenuation than massive ones, implying 
that the dust content is smaller, but the range of attenuation in 
more massive galaxies is larger. We also find that the dust at- 
tenuation in late-type galaxies increases as their NUV - r col- 
ors become redder, but early-type galaxies behave in the op- 
posite fashion. This probably indicates that the reddest early- 
type galaxies have little gas, ongoing SF and dust. 

Next we examine the total current SFR (averaged over the 
last 100 Myr). We find that GALEX is sensitive to very low 
SF levels, ~ lO~ 3 M yr -1 . Not surprisingly, galaxies with 
the lowest SFRs are primarily of early type. However, there 
is an interesting concentration of red dots at log SFR ~ 0.5. 
This "population" has lower concentration index and bluer 
NUV-r color than low-SF early-type galaxies. Inspection of 
the images of a subset of these galaxies with high metallicities 
reveals that about half of them are spirals (often disrupted), 
while most others show (low surface-brightness) disks. Some 
10% appear to be true ellipticals. This emphasizes the crude- 
ness of classification based on profile fitting. Finally, we 
note that we compared our SFRs with those derived from 
aperture-corrected nebular emission ( Brinchr nann et alJ2004l) 
and found excellent overall agreement. A more detailed com- 
parison of SF metrics will be presented elsewhere. 

The last row shows the ratio of current (last 100 Myr) to 
past-averaged SF (the "^-parameter"). Galaxies with constant 
SFR will have log b = 0, while those undergoing a burst can 
have log b > 0. We find that b is close to unity for low-mass 
galaxies but extends to 0.1 for m assive, late-types g alaxies, 
in agreement with the results of Brinchrna nn et alJ J2004I) . 
Early-types reach values as low as b rj 3 x 10~ 4 . The tight 
correlation between b and NUV - r color is most remarkable. 
It implies that the SF history of a galaxy can be constrained 
on the basis of the NUV - r color alone. 
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We have also examined the relation between SFR and FUV 
attenuation for the galaxies in our sample. We find that 
increasing SFRs are always associated with larger attenua- 
tions. At a fixed attenuation, however, early-type galaxies can 
have SFRs several orders of magnitude smaller than late-type 
galaxies. This suggests a different depend ence between dust 
content and SF in the two types of galaxies. iBuat et all (120041) 
also find that the dust attenuation increases with the amount 
of recent SF as estimated by the dust-corrected NUV flux. 

One might worry that because our sample has UV selection, 
the results shown in Figure |2] may not apply to the galaxy 
population as a whole. We have analyzed the relations be- 
tween the same galaxy properties as in Figure for the MIS 
fields, where the UV coverage of galaxies in the SDSS spec- 
troscopic sample is substantially more complete. We found 
exactly the same relations between SFR, b, A FUV and mass, 
concentration and color. The most notable difference is that 
the MIS data (by virtue of going deeper) contain a larger frac- 
tion of massive, concentrated red galaxies. Therefore, while 
the above results are useful for identifying various trends, they 
cannot be used to study the relative proportions of galaxies 
with different properties. 

Finally we perform a time-resolved analysis of starbursts. 
We select galaxies that have experienced starburst in which 
more than 5% of stellar mass (/-burst ^ 0.05) was formed over 
the last 100 Myr and 1 Gyr. Figure [3] gives the fraction of 
galaxies of a given mass that have undergone a recent star- 
burst. The dashed histogram corresponds to bursts within the 
last 100 Myr. For logM, < 9, the fraction is ~ 10%, but it 
declines to zero for the most massive galaxies. Over the last 
1 Gyr (solid histogram) some 20% of low-mass galaxies hav- 
ing experienced bursts. This fraction falls to 5% for larger 
masses. Qualitatively consistent results are obtained if we re- 
strict the analysis to the MI S fields. Similar results were ob- 
tained by Kauffman n et alJ i2003bl) . who analyzed starbursts 
over a 2 Gyr timescale. 

In this preliminary study of SF history using the UV pho- 
tometry from GALEX , we demonstrate the promise and point 
to possible limitations of this dataset when interpreted with 
galaxy models. As we gather more data we should be able 
to better explore the trends suggested here, and to get better 
statistics on low-mass galaxies. 
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